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Abstract 

I adapt a numerical method, previously applied to investigate the Yang-Mills vacuum wavefunctional, 
to the problem of extracting the effective Polyakov line action from SU(N) lattice gauge theories, with or 
without matter fields. The method can be used to find the variation of the effective Polyakov line action 
along any trajectory in field configuration space; this information is sufficient to determine the potential 
term in the action, and strongly constrains the possible form of the kinetic term. The technique is illustrated 
for both pure and gauge-Higgs SU(2) lattice gauge theory at finite temperature. A surprise, in the pure 
gauge theory, is that the potential of the corresponding Polyakov line action contains a non-analytic (yet 
center-symmetric) term proportional to |P| 3 , where P is the trace of the Polyakov line at a given point, in 
addition to the expected analytic terms proportional to even powers of P. 
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I. INTRODUCTION 



Consider a lattice gauge theory with gauge group SU(A0 on a periodic lattice of time extent 
N t , possibly containing matter fields and a chemical potential. If we integrate out all degrees of 
freedom under the constraint that Polyakov line holonomies are held fixed, then the resulting dis- 
tribution depends only on those Polyakov line holonomies or, more precisely, on their eigenvalues. 
The logarithm of this distribution is defined to be the effective Polyakov line action Sp. 

If the underlying lattice gauge theory in D = 4 dimensions has a sign problem due to a non-zero 
chemical potential, then Sp probably also has a sign problem. However, there are indications that 
the sign problem may be more tractable in Sp than in the underlying theory. Using strong-coupling 
and hopping parameter expansions, it is possible to actually carry out the integrations over gauge 
and matter fields mentioned above, to arrive at an action of the form 1 

3 

Sp = /3p£ £ [TiU}TiU x+ i + Trt/*Trt/; +i ] + k J>"Trt/, + e^Trl/}} , (1) 

x i=\ x 

where /3p, K are calculable constants depending on the gauge coupling, quark masses, and tem- 
perature T = l/Nt in the underlying theory. To minimize minus signs later on, the overall sign 
of Sp is defined such that the Boltzmann weight is proportional to exp[Sp], rather than exp[— Sp]. 
The Polyakov line holonomies U x E SU(A0 in (1) are also known as "effective spins." A path inte- 
gral based on an effective spin action of the form (1), for a wide range of /3p, K,fi, can be treated 
by a number of different methods, including the "flux representation" [2], reweighting [1], and 
stochastic quantization [3]. Even traditional mean field methods have had some degree of success 
in determining the phase diagram [4]. 

The problem, of course, is that strong lattice coupling and heavy quark masses lie outside 
the parameter range of phenomenological interest, and it is not obvious how to extract Sp for 
parameters inside the range of interest, even at ji = 0. There have been some efforts in this 
direction, notably the inverse Monte Carlo method of ref. [5], as well as early studies [6, 7] which 
employed microcanonical and Migdal-Kadanoff methods, respectively. There is also a strategy 
for determining the phase structure of lattice gauge theory from an effective spin theory, whose 
form is suggested by high-order strong-coupling and hopping parameter expansions [1]. Here, 
however, I will discuss a different approach to the problem, recently suggested in ref. [4], which 
will be illustrated for SU(2) pure gauge and gauge-Higgs theories. 



H. THE "RELATIVE WEIGHTS" APPROACH 

Let Sqcd be the lattice QCD action at temperature T = \/N t 'm lattice units, with lattice gauge 
coupling /3, and a set of quark masses denoted collectively m q . We set chemical potential /! = for 
now. It is convenient to impose a temporal gauge condition in which the timelike link variables are 
set to the unit matrix everywhere except on a single time slice, say at t = 0. In that case, Uq(x, 0) is 
the Polyakov line holonomy passing through the site (jc, t = 0). The effective Polyakov line action 



This is the action at leading order. For the effective action determined at higher orders in the combined strong- 
coupling and hopping parameter expansions, cf. [1]. 
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is defined in terms of the partition function 



Z(p,T,m q 



DU (x,0) J DU k DyDy e s ° CD 
DU (x,0) e Sp[Uo] , 



or equivalently 



exp 



Sp[U x ]}= / DU (x,0)DU k D\j7D\]f \ Y\8[U x -U (x,0)] \e s Q CD 



(2) 



(3) 



Because temporal gauge has a residual symmetry under time-independent gauge transformations, 
it follows that Sp[U x ] is invariant under U x — > g(x)U x g^(x), which means that Sp only depends on 
the eigenvalues of the Polyakov line holonomies. 

Now consider a finite set of M SU(A0 "effective spin" configurations in the three-dimensional 
cubic lattice V3 of volume L 3 , 



{{t/i ,all*eV 3 }, / = 1,2,..,m}. 



(4) 



Each member of the set can be used to specify the timelike links on the timeslice t = 0. Define 

& = J DU Q {x,0)DU k DwD W f^\Y[8[U^-U Q {xM\e SQCD , (5) 



i=l 



and consider the ratio 



exp 


S P [uU)] 


exp 





fDU (x,0)DU k Dxj7D ¥ {fL S[U^ - U (x,0)] } e s Qco 



fDUoix^DUkDWDxir \ U x S[ui k ) - U (x,0)} \e s oc° 



^ f DU (x,0)DU k D\}/Dy/ < 


[Ux^U^-Uoix^)} 


, e S QCD 


± f DUoix^DUkDxj/Dy < 


'Ux8[U^-U (x,0)][ 


, £ Sqcd 



(6) 



where in the second line we have merely divided both the numerator and denominator by a com- 
mon factor. However, by inserting this factor, both the numerator and denominator acquire a 
meaning in statistical mechanics, because the factor 5° can be interpreted as the partition function 
of a system in which the configuration of timelike link variables at t = is restricted to belong to 
the set {C/W, i = l,...,M}. This means that 



Prob[t/ W ] 



2? 



DUo(x,0)DU k Dxj7Dxj/ \ Y[8[U x j) - U (x,0)] \e s Q CD 



(7) 



is simply the probability, in this statistical system, for the 7-th configuration Uq{x, 0) = U^\x) to 
be found on the t = timeslice. This probability can be determined from a slightly modified Monte 
Carlo simulation of the original lattice action. The simulation proceeds by standard algorithms, 
for all degrees of freedom other than the timelike links at t = 0, which are held fixed. Periodically, 
on the t = timeslice, one member of the given set of timelike link configurations is selected by 
the Metropolis algorithm, and all timelike links on that timeslice are updated simultaneously. Let 
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Nj be the number of times that the ?'-th configuration is selected by the algorithm, and N tot = L/Af,-. 
Then Prob[*7 is given by 



Prob[£/ 

and this in turn gives us the relative weights 



exp 




exp 


~s P [u( k )}~ 



lim — — 

N, ot ^°° Ntot 



lim S 

AW-** 3 Nk 



(8) 



(9) 



for all elements of the set. A computation of this kind allows us to test any specific proposal for Sp, 
which may be motivated by some theoretical considerations. But it might also be possible, given 
data on the relative weights of a variety of different sets, to guess the action that would lead to these 
results. In this article we will consider sets of spatially constant Polyakov line configurations, 
and small plane wave perturbations around a constant background. This is already sufficient to 
determine the potential term in Sp, and to suggest the form of the full action. 

The method described above was proposed long ago [8] in connection with the Yang-Mills 
vacuum wavefunctional. Recently there have been some sophisticated suggestions for the form of 
this wavefunctional in 2+1 dimensions, and the technique was revived in order to test these ideas 
in ref. [9]. The main difference between the method as applied to vacuum wavefunctionals, and 
as applied to determining Sp, is that in the former case the simulation chooses from a fixed set of 
spacelike link configurations on the t = timeslice, while in the latter the choice is made from a 
set of timelike link configurations. 



A. Finite chemical potential 



Let Sq CD denote the QCD action with a chemical potential, which can be obtained from Sqcd 
by the following replacement of timelike links at t = 0: 



J QCD 



Sqcd\u (x,0) -> e^t/ (x,0),t/ t (x,0) -> e - N *U${x,0) 



(10) 



The corresponding Polyakov line action Sp is in principle obtained from (3), with Sq CD as the 
underlying action. Of course the integration indicated in (3) can so far only be carried out for 
strong couplings and large quark masses, but it is not hard to see that each contribution to Sp in the 
strong-coupling + hopping parameter expansion at /i = maps into a corresponding contribution 
to Sp by the replacement 



(11) 



It is reasonable then to suppose that this mapping holds in general, i.e. if we have by some means 
obtained Sp[U x , Ux] beyond the range of validity of the strong-coupling + hopping parameter ex- 
pansion, then the corresponding Sp is obtained by making the change of variables (11). There is, 
however, a possible source of ambiguity in this scheme (noted in [4]), coming from identities such 

as 



TrU 



t 



{TxU x ) 2 -TxUl 



(12) 
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in SU(3). One way around this ambiguity is to enlarge the range of Uq(x,0), allowing these 
variables to take on values 

U (x 1 0) = e i6 U(x) , (13) 

where U(x) is an element of SU (N). In other words, we allow the Uq(x, 0) links to take on values 
in the U (N) group, although it will be sufficient for our purposes to let 6 be Jt-independent. 2 
Suppose we are able to determine Sp for this enlarged domain of Polyakov line variables. Then 
Sp is obtained by analytic continuation, 6 — > — iN t fl. 

The essential point here is that if one can determine Sp by simulations of Sqcd at /i = 0, then 
this result can be used to determine Sp at finite chemical potential. If the sign problem is in fact 
tractable for Sp, as recent results seem to suggest, then this may be a useful way of attacking the 
sign problem in full QCD. 



B. Relative weights, and path-derivatives of Sp 

Let be the configuration space of effective spins {U x } on an L 3 lattice, and let the variable X 
parametrize some path {U X (X)} through . The method of relative weights is particularly useful 
in computing derivatives of the Polyakov line action 



dSp 



(14) 



dX 

along the path. To see this, we begin by taking the logarithm of both sides of eq. (9), and find 
Sp[U U) ]-S P [U {k) ] = lim {logA^-logA^ 



lim Jlog^L- log#-k (15) 



M f->°° 1 N tot N tot 



(From this point on we will drop the limit.) Now imagine parametrizing the effective spins by a 
parameter X; each value of X gives us a different configuration U X (X). Let the configuration 
correspond to A = Xq + AX, and U^ k > correspond to A = Xq — AX . Then 

ffggW) . i jiog^-iog^y (i6) 

V dX J x=Xo 2AX\ & N tol & N tot J 

However, rather than using only two configurations to compute the derivative, we can obtain a 
more accurate numerical estimate if we let X increase in increments of AX, e.g. 

X n = Xq+ (n - ^^—] AX , n=l,2,...,M , (17) 



2 It is also sufficient to restrict 9 to < 9 < 2n/N. The full range [0,27r] is redundant, because of the Zn center of 
SU(AT). 



5 



and use all of the M values obtained for N„ in the simulation. For AA small enough, the data for 
logN n /Ntot vs. X n will fit a straight line, and then we obtain the estimate 

( dS P [U x {X)] \ N n 
-p. « slope of log — vs. X n . (18) 

V dA /A=Ao N tot 
The procedure will be illustrated explicitly in the next section. 



HI. TESTING THE METHOD AT STRONG COUPLING 

The first step is to compute dSp/dX for a case where we know the answer analytically. As 
mentioned previously, Sp can be readily computed in the strong-coupling + hopping parameter 
expansion. We will consider here the case of pure SU(2) Yang-Mills theory at a strong coupling 
/3. If the lattice is N t lattice spacings in the time direction, then computing the diagrammatic 
contributions to Sp at leading and next-to-leading order in the strong-coupling/character expansion 
we find 



Sp 



X 1=1 



Nt 3 

E 

x i=\ 



(19) 



where 



P x = -TrU x 



j8p = 4 



N t 



(20) 



Let us first consider sets of spatially constant configurations with varying amplitudes in the 
neighborhood of P = Pq, i.e. 



(«) 



(P + a„ ) 1 + iyj 1 - (P + a n ) 2 a 3 
(n-^(M+l)^)Aa , n=l,2,...,M 



(21) 



so in this case a is the A parameter of the previous section. If we divide Sp into a kinetic and 
potential part, which in the case of (19) is 

Sp — Kp + Vp 



= ^fr E E ( p * p *+t - 2p x + p xPxi) 



x i=\ 
2 

x 1 

X 



(22) 



then dSp/da = dVp/dPo is giving us the derivative of the potential piece, which can then be 
reconstructed, up to an irrelevant constant, by integration. So the procedure for determining Vp 
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13=1.2, 12 J X 4 lattice 
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FIG. 1. The slope of the straight-line fit to the data shown gives an estimate for the derivative Lr^dSp/da 
of Sp with respect to the amplitude of spatially constant effective spin configurations. In this case, the 
derivative is evaluated at Pq = 0.5, for an underlying pure Yang-Mills theory at strong coupling value of 
J3 = 1.2, on a 12 3 x 4 lattice. 

(assuming it were not already known from the strong-coupling expansion) is to compute dVp/dPo 
numerically, fit the results to some appropriate polynomial in Pq, and then integrate the fit. 

Our sample simulation is carried out in pure SU(2) lattice gauge theory at coupling /3 = 1.2 
(well within the regime of strong couplings) on a 12 3 x 4 lattice with M = 20 sets of spatially 
constant configurations. Figure 1 shows the data for \og(N n /N tot ) plotted vs. (Po + a„)x spatial 
lattice volume (12 3 ), at Po = 0.5. It is clear that the data falls quite accurately on a straight line, 
and the slope gives an estimate for the derivative 

J_ ( dS P (U x (a)) \ 1 dVp(Po) 

Z 3 V da ) a=Q I? dP 

which can be compared to the value 6j5pPo obtained from the strong-coupling expansion. The 
derivative obtained from numerical simulation vs. Pq is plotted in Fig. 2, and it obviously fits a 
straight line. Therefore the potential Vp is quadratic in P x , and we find, at /3 = 1.2 

{0. 1721 (8) Ex \ p x relative weights method 
(24) 
0. 1710 J2jc \Px strong-coupling expansion 

where we have dropped, in the upper line, an irrelevant constant of integration. The small numeri- 
cal difference between the relative weights and strong-coupling results can probably be attributed 
to neglected higher order terms in the strong-coupling expansion. 3 

In order to investigate the kinetic term, we consider plane-wave deformations of spatially con- 

3 Statistical errors are estimated from best fit slopes obtained from eight independent runs. Where errorbars are not 
shown explicitly, in the two-dimensional plots shown below, they are smaller than the symbol size. 
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FIG. 2. A plot of the values for L^dSp/da vs. Pq. Each data point is extracted from a plot similar to the 
previous figure. Also shown are the corresponding strong-coupling values, and a best linear fit to the data 
points. 



stant configurations. The path through configuration space ^ is again parametrized by a, with 



U^=P^h + i^l-{P^) 2 a, 

Px = Po + a n cos(k-x) 
2% 

h = —mi , (25) 

where the {m,, i = 1,2,3} are integers, not all of which are zero. For this class of configurations 
we have, for the action (19) 

S P = p P L 3 ^3P 2 + ^£cos(*i) \ . (26) 

Since the deformation of the action is proportional to a 2 , it is natural to consider the derivative of 
Sp with respect to a 2 , i.e. 

^ = (27) 

and therefore we can choose to let a 2 , rather than a n , increase in equal increments, so that 
an = y/nAa. 

The numerical procedure is similar to the determination of the potential term: we compute the 
derivative L~ 3 dSp/d(a 2 ), at fixed Po an d k, from the slope of a plot of \og(N n /N tot ) vs. a 2 !?. 
Then these values for the derivative are plotted, at various values of Pq, against squared lattice 
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(3=1.2, 1 2 J X 4 lattice 



eg 

-5- 

tfi 
T3 

to 




1 1.5 2 2.5 3 3.5 4 



FIG. 3. Derivative of the action w.r.t. path parameter a 2 vs. squared lattice momentum. Data is taken at 
strong gauge coupling j3 = 1.2 for plane-wave deformations. Squares indicate the relative-weights values, 
while green dots are the values obtained from the strong-coupling expansion. 



momentum 

^ = 4£sin 2 (k-). (28) 
i=\ z 

The result, at Po = 0.5, is shown in Fig. 3, and we find, for a trajectory (25) at fixed k, 

1 dSp ~ -Ak 2 + B , (29) 



L 3 d(a 2 ) 
where 

A = 7.3(2) x 1(T 3 , B = 4.30(3) x 10~ 2 . (30) 

The simulation has also been carried out at other values of Pq, but the results are almost indistin- 
guishable from Fig. 3, and so are not displayed here. The important point, however, is that the 
path derivative (29) is Pq independent. 

Integrating with respect to a 2 , we find that along any path parametrized by a with fixed Pq 

S p [U x (a)]=L 3 {-Aa 2 k 2 L + Ba 2 +f(P )} , (31) 
where /(Pq) is a constant of integration, which can be determined from the data on the potential: 

f(P )=CP 2 , C = 0.0861 ±0.0004. (32) 
The next step is to express Sp along the path in terms of U x (or P x = jTril x ). From the definitions 



(25), (28), one easily finds that (31) can be expressed as 



Sp = 4A £ £ P x P x+i + \(B -6A)a 2 + (C- 12A)P^ 



JC ;'=1 



L 3 . (33) 



The constants B — 6A and C — \2A are, within statistical error, consistent with zero. So we will 
just drop these terms. Then along the trajectory the action has the form 

3 

S P = (.0292 ± .0008) £ £ P x P x+t (relative weights method) , (34) 

x i=\ 

and of course the natural conjecture is that this is the action itself, at any point in configuration 
space. Further checks would be to calculate numerical derivatives dSp / dX along other trajectories, 
to test the consistency of this conjecture. We don't really need to do that here, since the action 
at strong couplings is already known analytically, and is given in eq. (19) to leading and next-to- 
leading order in the strong-coupling expansion. At /3 = 1.2 we have, from eq. (19), that 

3 

Sp = .0285^^2 PxPx+i (strong-coupling expansion) , (35) 

x i=l 

which is a close match to what we have arrived at via the relative weights procedure. 

This is, perhaps, a lot of effort to derive a known result. We have gone through this exercise 
in order to illustrate the method, and to make sure, in a case where the answer is known, that the 
method actually works. 



IV. POTENTIAL V P IN PURE-GAUGE THEORY, WEAKER COUPLINGS 

We now reduce the lattice coupling of the underlying SU(2) pure-gauge theory, setting /3 = 2.2 
with inverse temperature N t = 4 in lattice units. At this coupling and temperature (which is still 
inside the confinement phase of the theory), the effective Polyakov line action Sp is not known. 

The easiest task is to determine the potential part of the action. For the purposes of this arti- 
cle, we define the kinetic part of the action to be the piece which vanishes for spatially constant 
configurations, while the potential part is local. With these definitions 

V P = ^y(U x ), (36) 

X 

and the function y(U x ) is determined by evaluating Sp on configurations U x = U which are con- 
stant in 3-space, i.e. 

r(£7) = ^S P (U) . (37) 

Then by definition the kinetic part of the action is 

Kp = Sp[U x ]-V P [U x ]. (38) 

In order to determine Vp, we consider as before the path through configuration space (21) 
parametrized by the variable a, and once again we can identify dSp/da with dVp(Po) /dPo as in 
(23). The derivatives are determined by the relative weight method described above, the depen- 
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dence on Pq is fit to a polynomial, and Vp is then determined, up to an irrelevant constant, by 
integration over Pq. 

Because the Z2 center symmetry is unbroken at /3 = 2.2 and N t = 4, and y(U x ) is a class 
function, it is natural to assume that Y(U) is well represented by a few group characters %j{U) 
of zero N-ality (j = integer for SU(2)), and the potential is analytic in P x . Surprisingly, this is not 
what is found. 

Figure 4(a) shows the data for the derivative 



D(P ) 



1 dV P 
~D~dP~Q 

1 dS P 
L 3 da 



(39) 



at /3 = 2.2 on a 12 3 x 4 volume, which, as in the strong-coupling case, extrapolates linearly to zero 
at Pq = 0. Also shown is a best fit of D(P) to the polynomial 



f(P) 



ciP + c 2 P 2 + c 3 P 3 



(40) 



with the best fit constants shown in Table I. What is initially a little troubling about this fit is that 
upon integration, and up to an irrelevant integration constant, we must have 



nPx) = \c\Pl + \ciPi + \c 3 pz , 



(41) 



which appears to violate center symmetry, i.e. y(P x ) = "V[—Px) for SU(2) gauge theory. Because 
of center symmetry, the character expansion of Y(P X ) contains only characters Xj with j = integer. 
It is a property of the SU(2) group characters that each %> can be expressed as a polynomial of 
order 2 j in P, containing only even powers of P for j = integer, and only odd powers for j = half- 
integer. Then if the character expansion of '1 / {P X ) is truncated at some j — jmax, the P-derivative 
is a polynomial in odd powers of P up to p 2 J ma ^ 1 . 

One might expect that i / (P x ) can be accurately approximated by a handful of group characters. 
However, the attempt to fit the data with only a few odd powers of P is unsuccessful, in the sense 
that each of the three fitting functions 



f(P) 



ClP + C 3 P 3 

ClP + C 3 P 3 +C 5 P 5 

ClP + C 3 P 3 + c 5 P 5 +c 1 P 1 



(42) 



corresponding to truncated character expansions with j max = 2,3,4, respectively, gives an un- 
acceptable fit, as seen in Fig. 4(b). The reduced % 2 values in the three cases are 440, 100,25, 
respectively. This is to be compared to the reduced % 2 = 3.2 for the fitting function (40). 



Potential fit 


c\ 




c 3 


4.61(2) 


-4.51(10) 


1.77(8) 



TABLE I. The constants ci_3 derived from a best fit of c\P + C2P 2 + C3P 3 to the potential data. 

All this seems to imply that y(P x ) has a term violating center symmetry, but of course that 
cannot be the case. In order that i / {P x ) is an even function of P x , it must be that the derivative is 
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(a) jj =gj? vs. P with a P 2 term in the fitting function. 
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(b) jj vs - P- Successive approximations without a P z term in 
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(c) test if ^p- is an odd function. 
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FIG. 4. Derivatives of the potential. Subfigure (a) shows the best fit to the data by a polynomial 
aP + bP 2 + cP 3 , while subfigure (b) shows a best fit by polynomials with two, three, and four odd powers of 
P, which are forms that might be expected from unbroken center symmetry, (c) is a test of whether dVp/dP 
is an odd function of P. Data for the derivative at values of Po < are multiplied by -1, for comparison with 
the data at Po > 0. (d) same data (and fit) as in subfigure (a), plotted in a different way. 
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an odd function, D(Pq) = —D(—Pq), which in turn means that the coefficient of the quadratic term 
in (40) must change sign when Pq — > —Pq. This is easy to check; we simply repeat the calculation 
with Pq < in (21), with the result shown in 4(c). Here the squares are the data for D(Pq) at Pq > 0, 
while the circles are data for (— 1) x D(Pq) at Pq < 0. The fact that the corresponding data points 
at ±Pq lie on top of each other means that the derivative is an odd function, and the potential itself 
is an even function of P x , as it must be. The conclusion, which follows from the best fit, is that 
over the full range — 1 < P x < 1 the potential, up to an irrelevant constant, is given by 

YiPx) = \c i Pl+ 1 -c 2 \P x \ 3 + 1 -c 3 Pl . (43) 

This function is non-analytic, because of the absolute value, but still center symmetric, with the 
constants given in Table I. It should be emphasized again that this potential cannot be approxi- 
mated very well by a simple sum of j = 0, 1,2,3,4 SU(2) group characters. Of course, any class 
function (including \P X \ ) can be approximated by a sufficiently large number of group characters, 
just as a step function can be approximated by a truncated Fourier series. But keeping only a rela- 
tively small number of group characters introduces "wiggles" in the approximation to the potential 
(which are seen in Fig. 4(b)) much like the truncated Fourier series does for the step function. 

So far we have only looked at a pure gauge theory in the confined phase, but it is also possible 
to compute y(P x ) in the deconfined phase using the same methods. In comparing the potential in 
the confining and deconfining phases it is useful to display the data in a slightly different way, by 
plotting the derivative dVp/d(P 2 ) vs. P 2 , i.e. 

LU(P 2 ) L?2PqcIPq ' 1 ; 

When the data is plotted in this way, a curious feature does show up. First, consider the confined 
phase. The data for the above derivative in the confined phase, at the same coupling /3 = 2.2 
and lattice volume as before, is shown in Fig. 4(d). In this plot, the best fit shown in Fig. 4(a) 
transforms to 

g(P 2 ) = \(c l +c 2 Vp 2 + c 3 P 2 ) 1 (45) 

with the same constants ci_3 shown in Table I, and this function is also plotted in Fig. 4(d). Note 
that if the potential didn't have a cubic term, then we would have to omit the term proportional to 
VP 2 . But then the data should fit a straight line in Fig. 4(d), which it quite clearly does not. 

Now we display corresponding data in the deconfined phase. Figure 5 shows the result for 
the derivative (44) at /3 = 2.4, again on a 12 3 x 4 lattice, which is well past the deconfinement 
transition. Note the peculiar "dip" near Pq = 0. Because of this dip, the polynomial form (40) to 
the derivative, which translates to (45) for dVp/d(P 2 ), cannot fit the data over the full range. It 
is consistent with the data away from the dip, i.e. at Pq > 0.1, and the resulting fit to data in the 
interval [0.1, 1] is also shown in Fig. 5. The relationship of the dip in the derivative near Pq = to 
the deconfinement phenomenon is not obvious to the author. 

Finally, it is important to ask whether the potential shown in Fig. 4 is dependent on the spatial 
volume. In Fig. 6 we show the previous data for the derivative of the potential, obtained on a 
12 3 x 4 lattice, together with data for the same observable obtained on an 8 3 x 4 lattice. It can be 
seen that the volume dependence is negligible in this case. 
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13=2.4, 12 J X 4 lattice 
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FIG. 5. Derivative of the potential in the deconfined phase. Note the dip in the data in the interval 
< F% < 0.1. The fit is to data at P n 2 > 0.1. 
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FIG. 6. A test of volume dependence of the potential at j3 = 2.2. Data for the potential derivative is displayed 
for lattice volumes 8 3 x 4 (open squares) and 12 3 x 4 (green circles). 



V. POTENTIAL V P IN SU(2) GAUGE-HIGGS THEORY 



We now add a matter field to the gauge theory, to see how this will affect the potential. To keep 
the computation requirements very modest, we consider a scalar matter field, in the fundamental 
representation, with a fixed modulus (i.e. a "gauge-Higgs" theory). For the SU(2) gauge group, 
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SU(2) Gauge-Higgs, p=2.2, 16 Lattice 



LU 




FIG. 7. Plaquette energy vs. gauge-Higgs coupling y at fixed j8 = 2.2, for the SU(2) gauge-Higgs theory 
with fixed Higgs modulus, showing a sharp crossover at y« 0.84. 



the matter field can be mapped onto SU(2) group elements, and the action can be expressed as 

S = P E lTr[UUUW] + yY,ln$Hx)U^)Hx+$)} ■ (46) 

plaq x,H 

There have been many numerical studies of this action, following the work of Fradkin and Shenker 
[10], itself based on a theorem by Osterwalder and Seiler [11], which showed that the Higgs region 
and the "confinement-like" regions of the /3 — y phase diagram are continuously connected. Subse- 
quent Monte Carlo studies found that there is only a single phase at zero temperature (there might 
have been a separate Coulomb phase), although there is a line of first-order transitions between the 
confinement-like and Higgs regions, which eventually turns into a line of sharp crossover around 
/3 = 2.775, y= 0.705, cf. [12] and references therein. At /3 =2. 2 the crossover occurs at y~ 0.84, 
as seen in the plaquette energy data shown in Fig. 7. There is also a steep rise in the Polyakov line 
expectation value as y increases past this point. 

Fig. 8(a) shows the potential derivative L~ 3 dVp/dPo vs Po, along with a best fit to the data, at 
/3 = 2.2 and y = 0.75, which is somewhat below the crossover, in the "confinement-like" regime. 
We compute this derivative, again in a 12 3 x 4 lattice volume, at both positive and negative values 
of Pq, to test for the presence of a small center- symmetry breaking term in the potential (which is 
not obvious in Fig. 8(a)). The data over the full range is fit to the form 

f(P) = c ; + c\P + c' 2 sign(P)P 2 + c ; 3 P 3 (47) 

which translates, upon integration, into a potential 

y(P x )=c' P x + y i Pl+ 1 -c' 2 \P x \ 3 + 1 -c' 3 Pl . (48) 

with a center symmetry breaking term c' P x . The constants obtained from the fit are shown in Table 
II. 

The slight asymmetry which breaks f(P) = —f(—P), and therefore center symmetry, is more 
evident when we expand the plot in the immediate region of Pq = 0, as in Fig. 8(b). It can be 
seen that the best fit through the data points does not go through f(Po) = at Pq = 0, but rather 



15 



Potential fit: gauge-Higgs model 






4 




0.025(1) 


4.70(2) 


-4.70(8) 


1.91 (7) 



TABLE II. The constants c' _ 3 derived from a best fit of c' + c\P + c' 2 sign (P)P 2 + c' 3 P 3 to the potential data 
of the SU(2) gauge-Higgs model. 



crosses the y-axis at a positive value f(0) = c' = 0.025. The line shown in Fig. 8(b) is taken from 
a best fit to the full range of data, not just the near Po = data. Since the underlying gauge-Higgs 
theory breaks center symmetry explicitly, a term linear in P x is of course expected. The coefficient 
Co = 0.025 of the symmetry breaking term is quite small, but the expectation value of the Polyakov 
line at 7= 0.75 is also quite small: (P x ) = 0.03 at these couplings and lattice size. 



P=2.2,i'=0.75, 12 J X 4 lattice 



2 
1.5 
1 

0.5 


-0.5 
-1 

-1.5 
-2 



data ■ 
fit 


w' 

/ 


7 

v 





-1 -0.8 -0.6 -0.4 -0.2 0.2 0.4 0.6 0.8 1 



(a) jj c Jj£ vs. P for the gauge-Higgs theory. 



P=2-2,Y=0.75, 12 J X 4 lattice 
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(b) A closeup near Pq = 0. 



FIG. 8. Derivative of the Polyakov line potential, per unit volume, with respect to the Polyakov line value 
P, for the SU(2) gauge-Higgs theory on a 12 3 x 4 lattice. Data is taken at gauge coupling j8 = 2.2 and 
gauge-Higgs coupling y = 0.75. (a) the data over the range — 1 < P < 1, together with the best fit; (b) the 
data in the vicinity of P = 0, also showing the fit in this region derived from the full range of data (i.e. same 
curve as in (a)). Note that the line through the data does not pass through the origin, which implies a small 
breaking of center symmetry. 



VI. PLANE-WAVE DEFORMATIONS 



We now return to the pure gauge theory at /3 = 2.2. So far the potential term Vp of the effective 
Polyakov line action has been determined, but the ultimate interest is in the full action. It was not 
very hard to extract this action from the log[N n /N tot ] data at strong couplings. Unfortunately it is 
not as easy to jump from the path derivatives to the full action at weaker couplings, simply because 
Sp is not so simple (and is not known in advance!). Nevertheless, knowledge of the action along a 
particular trajectory in configuration space does provide some information about the full action. 

As in the strong coupling case, we choose to investigate the derivatives of Sp along paths of the 
form (25), i.e. plane waves of fixed wavenumber and varying amplitude on a constant background. 
The method is the same as outlined in section III, but the result is different. At /3 = 1.2, it was 
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found that dSp/d(a 2 ) was linear in k 2 , and independent of Pq. That is not the case at /3 = 2.2. 
What happens in this case is shown in Fig. 9, where we display L~ 3 dSp/d(a 2 ) plotted against the 
magnitude of lattice momentum Icl = (k 2 ) 1 / 2 at fixed values of Pq — 0.1 and Pq = 0.8. It can be 
seen that the /^-dependence of the data in Fig. 9(a), at Pq — 0.1, is consistent with linear, while 
the ^-dependence in Fig. 9(b), at Pq = 0.8, seems to be quadratic. This can be seen from fits to 
a — bki in the former case, and to a — bk 2 L in the latter. This suggests a possible interpolating form 

1 dS P 



I?d( 



a. , „ D 2 



T -f(PQ)+cJk 2 +gP 2 , (49) 



whose ^-dependence would vary continuously from linear, as Pq — > 0, to quadratic, for k 2 L <C gP^. 
Fig. 10 is the same plot as Fig. 9(a), except that data obtained on both an 8 3 x 4 lattice and a 
12 3 x 4 volume are displayed together, and both sets of data points appear to have the same ki 
dependence. This is, of course, evidence of the insensitivity of our results to the spatial volume. 
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P = 0.8, B=2.2, 12 J X 4 lattice 



data 
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(b) Pq = 0.8 
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FIG. 9. Derivative of the action along a path of plane wave deformations, (a) Data at Po = 0.1 is consistent 
with a linear variation of the derivative with deformation lattice momentum ki, (b) data at Pq = 0.8 is 
consistent with a quadratic variation w.r.t. k^. 



If (49) is correct, then it ought to be consistent with the potential (43). This means that /(Pq) 
can be, at most, quadratic in Pq, so let us write 



1 dS P 

ZJdiV) 



bQ+b l PQ + b 2 PQ 2 +cJk 2 +gP 2 . (50) 



\a=0 



The constants shown are subject to three constraints by the potential, so if we insist on the potential 
(43) there are really only two independent constants. In order to derive those constraints, consider 
a very large lattice volume L 3 , such that k 2 can be made very small compared to gPy, but still 
non-zero, and we assume that {c\ , C2, C3} do not vary much with L (we have already seen evidence 
of this fact in Fig. 6). Then the kinetic term is negligible compared to the potential term, and along 
the trajectory (25), taking account of the spatial average (cos 2 k-x) av = \, we have 

1 dS P 1 1 „ 3 o 

- a ci + -c 2 Pq + -ciPq 2 . (51) 



L?d{a 2 )< 4 1 '2 ZUI 4 
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P = 0.1,13=2.2 



C/3 




FIG. 10. A check of insensitivity to lattice volume. Parameters are the same as in Fig. 9(a), but this time 
including data obtained on an 8 3 x 4 lattice volume (L = 8), in addition to data on a 12 3 x 4 volume (L = 12). 



Comparison with (50) in the k\ <C gP^ limit calls for identifying 

1 3 1 

bo=-c\ , b 2 = -c 3 , bi+cy/g= -ci ■ (52) 

Figure 1 1 show a best fit of the data to the form (50), with the best fit constants given in Table 
III. This is hardly a perfect fit through the data points, given the value of the reduced x 2 ~ 30. Still, 
except at very low £^,Pq, the fitting function gives a reasonable account of the dependence of the 
data on k\ and Pq. Table IV is a test of constraints, listing three combinations of constants which, 
according to the identities (52), should vanish. It is seen that the second and third combinations in 
the table are consistent with zero, and the first combination is very nearly so. 4 



Surface fit 


bo 


bi 


b 2 


c 


8 


1.105(14) 


0.85(17) 


1.365(56) 


-0.529(13) 


33(3) 



TABLE III. Fitting constants &o-2,c,g obtained from a best fit to the data points shown in Fig. 11, by a 
surface of the form (50). 



Constraints 






bi - |c 3 


-0.05(2) 


0.06(23) 


0.04(8) 



TABLE IV. The constraints (52) imply that the combination of constants in the second line of the table 
should vanish within errorbars, and the last line shows the actual values of these combinations, for the 
constants given in Tables I and III. 



4 All fits, and error estimates on fitting constants, are made using the GNUPLOT software. 
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(a) 




(b) 

FIG. 11. Two views, at different viewing angles, of the data (red crosses) for L dS/d(a 2 ) vs. lattice 
momentum ki and Polyakov line Pq, and the best fit (green surface) of the form (50) to the data. 
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VII. TOWARDS THE FULL ACTION 



The interesting question, of course, is what is the full action which gives rise to the variation 
(50) along the path, with the given potential (43). We begin by noting that, with the constants 
shown in Tables I and III, the action 



x I 1 



S P = 2c ^£P x Q xy P y - £ y^P, 2 1 + E(^iP 2 + l -c 2 \Pl\ + l -c,P A x 

= K P + Y d nPx) (53) 



X 

where Kp is the kinetic term 



K P = 2c I "£P x Q xy P y - £ <Jg~P 2 P 2 1 (54) 

[xy X J 

and 

Qxy=(VR) xy 
R X y=(-^l)xy+gP^xy 

= £(25^ - 8 x , y+i - 8 x+i ) +gPl8 xy , (55) 

;=1 

gives the known results for the potential (43) and for the variation of Sp with a 2 (50) along the 
paths of plane wave deformations (25). The operator is the usual lattice Laplacian operator, 
and Q has the spectral representation 

e=L(v / *£+^) !*><*! 

0^ = ^1,(^ + 8^) e iHX ~ y) > (56) 

where Y.k 1S shorthand for the sum over lattice wave vectors with components kj = (2^/L)m,-, 
and lattice momentum ki has been defined previously in (28). The ket vectors \k) correspond to 
normalized L~ 3 / 2 exp[ik- x] plane wave states. 

For the paths (25), set P x = Pq + acos(k ■ x), and compute the resulting action on such config- 
urations up to leading order in a 2 . Using the spectral representation for the operator Q, a short 
calculation gives, up to 0(a 2 ), 

S P = L 3 r(P ) +a 2 L 3 jici + C-c 2 -c^)P + ^c 3 P$ + c^Jk 2 +gP 2 ^ . (57) 

Applying the identities (52), which are reasonably well satisfied by the data, this becomes 

S P = L 3 r(P )+a 2 L 3 |Z7 + Vo + b 2 Po +c^/k 2 +gP 2 ^ . (58) 

So we find that for constant configurations (a = 0), the action is simply the known potential, i.e. 
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Sp = L 3 i / (Po), while the path derivative is 

1 dS P 



L 3 d(a 2 ) 

in complete agreement with (50). 



b + b l P l +b 2 P£+c x /ki + gPj , (59) 



Denote by P av and AP 2 the lattice average value and mean square deviation, respectively, of a 
given Polyakov line configuration. It is clear that for the paths (25) considered so far, Pq = P av . 
One further generalization, which will not affect agreement with the data so far, is to allow the 
kinetic term to also depend on AP 2 , i.e. 5 

K P = 2chp x U-V 2 + gP 2 , + g'AP 2 ) P y -^^gP 2 , + g'AP 2 pA (60) 

[ xy W x J 

It is not hard to see that the 0(a 2 ) contribution that would arise from the a 2 -dependence of the 
square root terms also selects, at this order, the constant a 2 -independent part of P x and P y . In that 
case ki = 0, and this contribution to the 0(a 2 ) part of the kinetic term vanishes. 

In order to investigate the possibility of a AP 2 -dependence a little further, let us consider tra- 
jectories consisting of plane waves, of varying amplitude A, with P av = 0, i.e. 

P x =Acos(k-x) , (61) 

and study the derivative L~ 3 dS p /dA evaluated at A = Aq. To compute this derivative by the relative 
weights approach, we construct a set of configurations 



Px = A n cos(k-x) 



A n =A Q +(n-^{M+\)^AA , n=l,2,...,M 



In 

T 

and proceed as before. The conjectured action is 



h = —mi . (62) 



Sp = 2c I ^P x U-V 2 +gP 2 , + g'AP 2 ) P y -£^gP 2 v + g>AP 2 P } 

[ xy 7 *y x 

+ H(\^ P x+lc2\P^\ + \c 3 P } 



(63) 



5 



A generalization of (53) which does not work is the replacement of Pq by P x in (53) and (55). This leads to additional 
contributions to dSp/d(a 2 ) which spoil the agreement with (59). 
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FIG. 12. Variation of Polyakov line action with Polyakov line amplitude, L~ 3 dSp/dA evaluated at A = Aq, 
for Polyakov line configurations proportional to plane waves P x = A cos (k • jc) , as a function of Ao and lattice 
momentum k^. Red crosses are data points, and the green surface is a best fit to the data by the analytic 
form (64). 



whose path derivative is 
1 dS P 



L 3 dA \ a=Aq 2 

1 



l - Cl A Q + A2Ac 2 Al + .315c 3 A 3 + 2cA (^J k 2 + ^g'A 2 - \J\g^ 



+ -cg'Al I - ' - -== | (64) 

2 \ iki + yAi Jy^ t 



Taking c and ci_3 as given in Tables I and IV, there is only one free constant left to fit the 
data, and the best fit, shown in Fig. 12, is obtained at g' = 3.45(4). Once again, this plot should 
not be interpreted as a perfect fit through the data points within errorbars, given that reduced 
X 2 ~ 45. On the other hand, with only one fitting constant, the expression (64) does seem to give a 
quite reasonable account of the dependence of the data on Ao and Ul, despite the highly non-local 
expression AP 2 introduced into the kinetic term. 



6 The numbers multiplying c\ , C2 , C3 are the lattice averages of cos 2 (k ■ x), \ cos 3 (Jfc • x) \ , cos 4 (k ■ x) respectively. These 
numbers are almost independent of the wavenumber k on finite lattices, so long as k 7^ 0, and converge rapidly to 
the infinite volume limit as lattice volume increases. 
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VIII. CONCLUSIONS 



I have presented a method for computing derivatives dSp/dX of the effective Polyakov line 
action along any given path through field configuration space, parametrized by the variable X. 
The technique is easily implemented in a lattice Monte Carlo code by simply replacing updates of 
timelike links, on a single timeslice, by a Metropolis step which updates that set of links simulta- 
neously, and the potential part Vp of the effective Polyakov line action can be readily determined, 
for any given lattice coupling, temperature, and set of matter fields, up to an irrelevant constant. It 
is also possible to determine, from the derivatives, the action Sp along any given trajectory in field 
configuration space. 

The method has been applied here to SU(2) lattice gauge theory, both without and with a scalar 
matter field. At a strong coupling (/3 = 1.2) and finite temperature, the method easily determines 
the effective Polyakov line action, which we have checked against the known result derived from 
a strong-coupling expansion. At a weaker coupling (/3 = 2.2 on a 12 3 x 4 lattice), where the 
Polyakov line action is not known, it has been shown that, up to a constant, the potential term has 
the form 



with coefficients given in Table I. The center- symmetric but non-analytic cubic term comes as a 
surprise; to the best of my knowledge such a term has not been anticipated in previous studies. It 
would be interesting to study the evolution of the above potential as /3 and N t vary. Addition of 
a scalar matter field in the underlying lattice gauge theory introduces a center symmetry breaking 
term into the potential which is linear in P x , with a coefficient reported in section V. 

Data has also been obtained from small plane-wave deformations around a constant Polyakov 
line background (Section VI), and for Polyakov lines proportional to a plane waves with variable 
amplitude (Section VII). It was found that the action (63) is consistent with the results that have 
been found so far, and at this point we may conjecture that (63) approximates the desired full 
Polyakov line action. Of course, the kinetic term in Sp could easily have a more complicated 
form than what is suggested in (63), and therefore this conjecture needs to be tested on more 
complicated, non-plane wave configurations. Those tests, and the extension to the SU(3) group, 
would be the obvious next steps in the approach introduced here. 
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